Thermally driven ballistic rectifier 
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The response of electric devices to an applied thermal gradient has, so far, been studied almost 
exclusively in two-terminal devices. Here we present measurements of the response to a thermal 
bias of a four-terminal, quasi-ballistic junction with a central scattering site. We find a novel 
transverse thermovoltage measured across isothermal contacts. Using a multi-terminal scattering 
model extended to the weakly non-linear voltage regime, we show that the device's response to a 
thermal bias can be predicted from its nonlinear response to an electric bias. Our approach forms 
a foundation for the discovery and understanding of advanced, nonlocal, thermoelectric phenomena 
that in the future may lead to novel thermoelectric device concepts. 

PACS numbers: 72.20.Pa, 73. 23. Ad, 73.40.Ei 
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Nanoelectronic devices with dimensions smaller than 
the characteristic scattering lengths for electrons are eas- 
ily driven into the nonlinear regime, and show electron 
dynamics that depend on device symmetry. For instance, 
clear rectification effects have been observed in asym- 
metric microjunctionai""— , quantum dots^ and quantum 
point contact systems^i^. These systems exhibit non- 
equilibrium effects that occur in the nonlinear regime of 
ballistic transport. 

At the same time, there has been considerable progress 
in the use of mesoscopic conductors for thermoelec- 
tric devices such as nanoscale coolers and ultrasensitive 
thermometers^. At a more fundamental level, studies 
of heat flow in phase-coherent systems allow one to in- 
vestigate basic properties of quantum transport when 
electric and thermal gradients act jointly on the same 
systemp^.. Specifically, thermoelectric effects are observed 
when a voltage drop builds up in response to a temper- 
ature difference^iiiS. Thus far, most of the literature has 
investigated thermopowers in two-terminal arrangements 
(for exceptions see Refsi^iii). 

Additionally, for the past four decades there has been 
sporadic research on so-called anisotropic thermoelec- 
tric materials. These are materials that exhibit off- 
diagonal elements in their electrical conductivity, ther- 
mal conductivity, or thermopower. Typically, anisotropic 
thermoelectrics comprise of either stacked multilayers 
of alternating materialai^ii^, or anisotropic crystalline 
materialsi^. Such systems, however, can only change 
the anisotropic properties of a given system by rotating 
the stacked multilayer or crystal orientation. This effec- 
tively limits the range of anisotropic properties that can 
be studied. 

Here, we address transverse thermoelectric effects in 
a quasi-ballistic device fabricated in a two-dimensional 
electron gas (2DEG), Fig. [TJa), which offers the great 
advantage that the scattering properties of the system 
can be controlled during fabrication. Specifically, we use 
a four-terminal microjunction with a central, asymmet- 



ric antidot to scatter electrons, Fig. [TJb), and perform 
thermoelectric measurements at a cryostat temperature 
of 0.74 K. In a similar device, previous observations have 
shown that a lateral AC voltage leads to a rectified, trans- 
verse response^. Here, we instead apply a lateral temper- 
ature difference across the junction through terminals 1 
and 2, and observe a nonlocal thermovoltage between ter- 
minals 3 and 4, both of which are unheated. We note that 
similar thermovoltages have been detected previously in 
multi-terminal 2DEGs subjected to magnetic fields nor- 
mal to the system (the Nernst effectiiii^'i^) . In contrast 
to these works, the effect reported here is purely thermo- 
electric in nature and only due to the symmetry breaking 
property of the central scatterer. Furthermore, we es- 
tablish a theoretical model based on the multi-terminal 
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FIG. 1. (a) Schematic of the wafer structure with a 2DEG 
located in the GalnAs quantum well, (b) Schematic and SEM 
image of a four terminal antidot junction similar to the one 
reported on here. The four circles indicate connections to dis- 
tant Ohmic contacts, not shown. The circuit diagram on the 
right shows the heating circuit used to heat the side channels; 
the configuration shown here is used to heat the right side 
channel. The lateral and transverse voltage drops are given 
by Vi2 and V34, respectively. The actual four-terminal junc- 
tion measured here had terminal widths of di « 500 nm and 
^2 ~ 450 nm, and an antidot with a base and height of about 
550 nm and 380 nm, respectively. The side heating channels 
are 5 /^m wide, 30 /^m long, and ~ 4 /im from the scatterer. 
Note that the side channels are not to scale. 
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FIG. 2. (Main figure) First four iiarmonic transverse voitage 
responses to a temperature increase in terminal 1. Tlie solid 
red line is a parabolic fit to the 2io curve at low voltages, and 
is given by V'3'4"' = (1.9V"^)Fh . (Top inset) Schematic of the 
measurement configuration. (Bottom inset) Measured lateral 
voltage drop across the junction during heating. 



scattering approachiii^ extended to the weakly nonlin- 
ear regime (up to quadratic order in voltage and linear 
order in temperature shifts). We then show that this 
model can predict the transverse thermoelectric response 
to a lateral thermal bias from nonlinear voltage measure- 
ments similar to those reported in Ref.^. 

Our device consists of an InP/Gao.23lno.77As 2DEG 
wafer, see Fig.[TJa), patterned using the same techniques 
as in Ref.l^. The antidot junction itself and its dimen- 
sions are shown in Fig. [Ijb) . Hall and Shubnikov de Haas 
measurements give a carrier concentration and mobility 
of 2.93 X 10^^ cm~^ and 1.16 x 10^ cm^/Vs respectively. 
The resulting mean free path, 3.3 ^m, exceeds the char- 
acteristic length of the antidot region, 0.5 pxa.. Therefore, 
electrons within the antidot junction are in the ballistic 
regime. 

For our thermoelectric measurements, we are primar- 
ily interested in the transverse voltage response, V34 = 
^3 ~ ^j due to a change in the side terminal tempera- 
ture(s), of either 0i, Qi, or both. Here Vi and Qi are the 
shifts of the electrochemical potential and temperature in 
the i terminal away from the common electrochemical 
potential, /.i/e, and background cryostat temperature, d. 
To increase either B\ or ^2, we apply two 37 Hz heating 
currents, which have a relative phase difference of 180°, 
to the channel to be heated. Figure [Ijb) shows the cir- 
cuit diagram for heating the right side channel. Using 
P\ as a voltage probe, the relative amplitudes of vj^ and 
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are tuned so that the resulting la; voltage at terminal 
1 was minimized. Any remaining DC offset could have 
been minimized as well using a DC source to shift both 
heating voltages; however, no offset was needed. This 
technique allows us to apply a thermal bias without an 
electrical bias. 
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FIG. 3. (Main figure) l'^* and 2"^* harmonic transverse re- 
sponses, V34, to a laterally applied 37 Hz voltage 1^12. Also 
shown is a parabolic fit to the 2lo curve at low voltages, (inset) 
Circuit schematic for the applied voltage measurement. 



In Fig. [2] we show results from an experiment where a 
temperature shift d\ was applied to terminal 1, and the 
first four harmonics, with respect to the heating volt- 
age, of V34 were measured. Since the thermal gradient 
is generated from Joule heat, which is primarily pro- 
portional to the square of the applied heating voltage, 
the dominant 2aj response implies that the transverse re- 
sponse is proportional to the applied temperature gra- 
dient, V34 oc l/J oc 01. We also measured the first and 
second harmonics of the lateral voltage drop V12 using 
the probes labeled P\ and P2 (see inset to Fig. [5]). The 
lateral voltage by itself, V\2 < 4 /itV, is far too small 
to cause the observed transverse voltage, as can be seen 
from the device's measured Iw and 2uj responses to an 
applied 37 Hz voltage, see Fig. [31 

At low voltages, V34 in Fig.[2]is linear in 9i and can thus 
be interpreted as a thermovoltage. Remarkably, however, 
this thermovoltage is across two unhealed contacts that, 
based on device geometry, can be assumed to be isother- 
mal. We therefore propose that ^34 is due to the differ- 
ence between two or more thermovoltages in the junction. 
Sources of thermovoltages that are easily identifiable are 
the four quantum point contacts (QPC) defined by the 
antidot and the surrounding channel walls. In the general 
case of zero symmetry, each QPC thermovoltage would 
be unique, resulting in all four terminals being at dif- 
ferent potentials. We will return to this interpretation 
below. 

Given the intended left-right geometric symmetry of 
the junction one may intuitively expect that ^34(^1, 62) — 
V34{92,0i). When the lateral temperature differential is 
reversed, however, measurements show a sign reversal 
in V34 (Fig. U]). Furthermore, when both terminals are 
heated simultaneously, the response tends toward zero. 
We will return to these unexpected observations shortly. 

We now aim to explain our results using the multi- 
terminal thermoelectric scattering approacbiii^ ex- 
tended to the weakly nonlinear transport regime^S. We 



consider a four-terminal junction with a central scatterer, 
as in Fig. [IJb) and Ref J^, but without assuming any par- 
ticular model for the scattering matrix. Based on the lin- 
ear temperature and quadratic voltage responses shown 
in the experiment (see Figs. [2] and [3]), we express the cur- 
rent la through lead a up to linear order in temperature 
and quadratic order in voltage differences, Vap = Va — Vp: 
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[GapiVap) + Fapiy^pf + LapiOa - Op)] . (1) 



Importantly, Eq. ([T]) is manifestly gauge invariant, a 
property recently emphasized for mesoscopic rectifiers^i. 
In the low temperature limit, the linear transport co- 
efficients read, 
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This approximation is valid provided the energy variation 
of the transmission probability, tap , between leads a and 
/3 is weak and the Fermi energy, ep, is larger than the 
background temperature, 9. This amounts to neglecting 
band edge effects. Hence, the transmission and its energy 
derivative, t'^p, are evaluated at ep- 

The leading order rectification term Fap is a compli- 
cated function of t'^a and of the transmission deriva- 
tive dtap{E,{Vi})/dVa\cq evaluated at equilibriun>2£. 
Knowledge of the latter requires a self-consistent determi- 
nation of the screening potential, U{{Vi}), inside the con- 
ductor. However, within a qualitative discussion we can 
neglect this term since it is proportional to the character- 
istic potential, Ua — {dU/dVa)\ccr^, whose strength can 
be estimated within mean-field theory from u ~ C^/C, 
where C~^ — {e'^D)~^ + C~^ , C is a capacitance co- 
efficient that measures interactions, and D denotes the 
density of states associated with the scattering states. In 
widely open systems one has C ^ e^D (noninteracting 
limit), and the rectification term becomes 



FaP ^ -J^t'api'^F)- 



(4) 



To find V34 as a function of the lateral voltages {Vi and 
V2) and temperatures {0i and 62) we recall that terminals 
3 and 4 act as voltage probes. We then solve the two 
conditions I3 = I4 — with the aid of Eq. ([T]) in the cold 
isothermal case (^3 = ^4 = 0), that is, we assume that 
the channel thermal conductances are much greater than 
the thermal conductances across the antidot region. For 
definiteness, we neglect the products ViVi^4\ since in our 
experiment, V34 is much smaller than V^i2. Then, 
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FIG. 4. 2a; transverse response due to heating both side chan- 
nels simultaneously and individually. The vertical dotted lines 
indicate the approximate onsets of nonlinear regimes. 



where we have defined, 

. ^31^42 — ^41^32 
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Bj = 



(til + ^42)^3^ — (^31 + t32)t4j 
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(7) 
32)- 



3=1,2 



withZ) = {t3i+t32)it4l+t42) + t34{t4i+t42)+t43it3i+t 

Here ge — e/2 and gg = Tr'^kg6/3e. 

Interestingly, Eq. [S] predicts that the coefficients Bj 
that describe the non-linear voltage response also de- 
scribe the transverse response to a laterally applied tem- 
perature differential. We can utilize this relationship 
to extract the expected thermopower from the electrical 
measurements in Fig. |3l and the electron temperature 
rise from the thermal bias measurements in Fig. [2] 

From the parabolic fit in Fig. [31 and Eq. ([5]) , we de- 
termine Bi — 2.02 X 10^" J~^. Inserting this value 
back into Eq. ^ yields a transverse thermopower of 
S'34,12 — Bigg — 594 nV/K. We will address the mag- 
nitude of this thermopower in a moment. 

With 5*34^12 in hand, we can then estimate 9i as a 
function of Vh- Using the parabolic fit from Fig. [21 we 
find 61 w (3.2 X 10^ K/V^)V^2^ which is vahd for I4 < 
0.2 mV. At the upper limit of the fit, Vh = 0.2 mV, 
this yields 6*1 « 130 mK. These temperatures agree well 
with a simple diffusive-heat finite-element simulation of 
the electron temperature in the side channel: 6*^™' ~ 
(2.5 X 10^ K/V^)V^2 foj, y^ ^ 0.2 mV. 

To offer a more conceptual description of our observa- 
tions, we now interpret the effective thermopower of the 
junction as the difference of QPC thermopowers between 
neighboring terminals. For the configuration where heat 
is applied to the right channel, we interpret the effective 
thermopower as the difference between the top-right and 
bottom-right QPC thermopowers. Correspondingly, we 
interpret the thermopower for heating of the left channel 
as due to the left-hand pair of QPCs in the junction. To 



check whether this interpretation is reasonable, we com- 
pare the magnitude of 5'34^i2 to the difference between 
two ideal QPC thermopower peaks: 
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where A'^rT.s) ^ 1 represent the integer number of con- 
ductance modes in the respective top and bottom QPCs. 
Using the measured two-terminal conductance between 
terminals 1 and 2, G'^l"'"'- = 8.73 x lO^^^r^-i, we estimate 
the number of modes to be around 11. For mode num- 
bers between 10 and 13, Eq. ([5]) predicts two-QPC ther- 
mopowers with a magnitude between 0.4 and 1.5 /iV/K, 
which fits well with S'34_i2. Here it is important to note 
that the true two-QPC thermopower is a function that 
strongly oscillates between positive and negative values^, 
and Eq. 1^ only predicts the range of reasonable values. 
However, the oscillatory nature of this thermopower of- 
fers a simple explanation for the sign change in the ther- 
mopower observed in Fig. |4j because the sign of the ther- 
mopower depends sensitively on the number of occupied 
modes in each QPC, even a very small left-right asym- 
metry in the width (occupation number) of the QPCs in 
our ballistic junction can lead to a sign change. 

Our conceptual two-QPC model also offers an expla- 
nation for the nonlinear changes in slope of the ther- 
mopower as a function of heating power (see vertical, 
dashed lines in Fig. |4]), which have been seen in simi- 
lar two-QPC systems^. As the temperature drop across 
a QPC increases, thermal broadening about the Fermi 
energy accesses neighboring thermopower peaks and in- 
creases the QPCs thermopower. The downward change 
in slope in the heat-right curve can then be interpreted 
as an increase in the number of accessed thermopower 
peaks in the top right QPC as a function of 0i. Similarly 
for the heat-left configuration, the two nonlinear points 
indicate increases in the number of accessed thermopower 
peaks in the bottom left QPC. 

It is important to note that there are in fact four in- 
teracting QPCs in the junction, a point the two-QPC 
model partially neglects. To properly interpret the ob- 
served sign change in Fig. 2] then, we refer to ([S]), which 
simultaneously accounts for all four QPCs. Focusing only 



on the response to temperature changes, the 2a; response 
can be decomposed into components that are symmetric, 
and antisymmetric, 5'; 
i,c'2) -^ (^2,^i): 
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^, under the transformation 



1^34 = S: 
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where S^^ = Bi + B2 and S^^ = Bi - B2. When 
both terminals are heated simultaneously and equally, 
the main contribution stems from the symmetric compo- 
nent, which is generally nonzero. When 9i and 62 are not 
equal and exchanged, the symmetric component remains 
unchanged while the antisymmetric component reverses 
sign. Note that S^^ exists only for scatterers showing 
some kind of asymmetry in the transmission derivative 
coefficients (^3^ 7^ t'^2 or t'^i ^ ^42). Since we observe a 
complete change in sign, our device must have a strong 
asymmetric component. 

We have shown that a four-terminal junction with a 
symmetry-breaking scatterer can be used to generate 
a thermovoltage between two cold isothermal contacts. 
We have put forward a simple noninteracting model for 
multi-terminal thermoelectric transport that accounts for 
most observed effects. Given the great control over de- 
vice geometry and symmetry afforded by nanofabricated 
devices, our approach can be used to explore advanced 
thermoelectric functionalities — such as electronic ther- 
mal rectifiers or more efficient thermoelectric energy con- 
verters — due to the opportunity to separate heat and 
charge fiow in multi-terminal devices. 
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